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ABSTRACT 


The majority of the manned and unmanned missions in space for at 
least the next ten years will be composed mainly of earth orbiting 
satellites and space stations. Although there are many environmental 
aspects of space which must be explored and elaborated upon before the 
successful launching of such missions, few are so critical as the 
radiation environment which will be encountered. 

The radiation environment for earth orbiting vehicles will consist 
primarily of (1) Van Allen particles, (2) solar flare particles, and 
(3) galactic cosmic particles. Of these three, the Van Allen belt 
particles present the greatest hazard. 

This report summarizes the knowledge of the energy spectrum and 
distribution of the Van Allen proton belt as it is known today. A 
computer program has been developed for calculating Van Allen particles 
incident upon vehicles in earth orbits. Based upon the 1963 FLUX program 
of Me II wain, we find that a vehicle in a highly elliptic orbit may 
encounter as many as 1.1 x 10 8 protons /cm 2 / orbit with 40 Mev < F < 110 
Mev. In the magnetic equatorial plane, the flux of protons with E > .5 
Mev may be as high as 2.2 x 10 7 protons/cm 2 /sec. > 
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AN ENVIRONMENTAL MODEL FOR VAN ALLEN BELT PROTONS 

SUMMARY 


The majority of the manned and unmanned missions in space for at 
least the next ten years will be composed mainly of earth orbiting 
satellites and space stations. Although there are many environmental 
aspects of space which must be explored and elaborated upon before the 
successful launching of such missions , few are so critical as the 
radiation environment which will be encountered. 

The radiation environment for earth orbiting vehicles will consist 
primarily of (1) Van Allen particles, (2) solar flare particles, and 
(3) galactic cosmic particles. Of these three, the Van Allen belt 
particles present the greatest hazard. 

This report summarizes^ the knowledge of the energy spectrum and 
distribution of the Van Allen proton belt as it is known today. A 
computer program has been developed for calculating Van Allen particles 
incident upon vehicles in earth orbits. Based upon the 1963 FLUX program 
of Mcllwain, we find that a vehicle in a highly elliptic orbit may 
encounter as many as 1.1 x 10 s protons /cm 2 / orbit with 40 Mev < F > 110 
Mev. In the magnetic equatorial plane, the flux of protons with E > .5 
Mev may be as high as 2.2 x 10 7 protons/cm 2 /sec. 


I. INTRODUCTION 

Before 1958 the Van Allen radiation zones were not known to exist, 
and the only forms of radiation expected outside of the earth’s atmosphere 
were the low flux, high energy, cosmic rays and occasional solar flares. 
Although little information was available about solar flares, they were 
known to be composed of high energy particles and to correlate temporally 
with gigantic eruptions on the photospheric surface of the sun. 

The first American satellite, Explorer I, carried a Geiger counter 
for the purpose of obtaining further information on the flux of cosmic 
ray particles. The experiment worked well until the satellite reached 
about a thousand kilometers altitude, at which point the counting rate 
suddenly decreased. The experimenters (Van Allen, Mcllwain, Ludwig, and 
Ray) were suspicious of what was happening, but did not publicize their 



suspicions until after some Explorer III data were examined. When the 
same phenomenon was observed on Explorer III, Van Allen announced 
(Ref. 38) that the apparent decrease in the Geiger counter was due to 
saturation by a high flux of particles. 

These particles, which are trapped by the magnetic field of the 
earth, present a radiation problem which must be overcome before 
extended earth orbital missions of manned spacecraft can be performed. 

In addition, the environment to which unmanned satellites are exposed 
may decrease the lifetime of on-board experiments, power supply, telemetry, 
etc. Thus, extensive environmental studies of the radiation zones should 
be conducted before the development of sophisticated long-lifetime 
satellites . 

This report summarizes the known facts and suggests a possible 
environmental model for Van Allen belt protons. 


II. MOTION OF CHARGE PARTICLES IN THE MAGNETOSPHERE 

The motion of charged particles trapped in the earth 1 s vicinity by 
the geomagnetic field is influenced by this field. Perturbations caused 
by an electric field imposed by the earth’s equatorial ring current also 
affect the orbital motions of these charged particles. 

We may consider a charged particle moving in this magnetosphere^ as 
being subjected to the geomagnetic field, B, and a variable field, E, 
due to the high altitude (~6 R e ) ring current. (R e = geocentric earth radii). 
The force, 1?, exerted on the particle moving with velocity, v, is 


— > — » — > 

F = q (E + v X B) , 


where q is the total charge of the particle. We will here consider^ only 
the effects of the geomagnetic field on the particle. The force, F, is 
normal to the plane formed by v and B; thus, the particle is induced, 
assuming that no other forces are present, to move in a circle with no 
loss in kinetic velocity. 

If we consider the portion of v normal to the magnetic field as v^ , 
we may say that 


mVj 2 

— = qv. B, 

r * 


2 



where r is the Larmor radius, the radius of the circle in which the 
particle moves, or 


r 


mVj_ 

q~B 


The Larmor frequency, oj, is then 



The particle is now in space moving in a circle which is dependent upon 
the magnetic field strength, the particle's charge, mass, and velocity 
normal to the magnetic field. When the particle begins to orbit, the 
center of the orbit is located on a particular "line of force" of the 
magnetic field. This specific line of force then remains the center of 
gyration or "guiding center" of the particle, unless the particle is 
subjected to some external force. The orbit of particles trapped in the 
geomagnetic field is subjected to perturbations caused by the presence 
of the ring current, as well as inhomogeneities and slow variations of 
the field itself. 

The component of velocity, v M , which is parallel to the magnetic 
field must also be considered. This component of the particle's velocity 
causes the particle to oscillate between "mirror points" in the northern 
and southern hemispheres. The particles are reflected at these points 
of equal magnetic field strength. When the particle is at the magnetic 
equator of the earth, it is at the point of minimum field strength on 
that line of force. As it moves either north or south along that line 
of force, it encounters a stronger field which in turn exerts a force, 

F„, on the particle, which is opposed to v„. 

This force is -dependent upon the magnetic moment, |_i, of the particle. 


where 


or 


r 2 q oj 


H = 



B 
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Finally, the retarding force on a particle moving along a line of force 
into a region of more intense magnetic field is given by 

F u = - p (V B), 

where the (V B) denotes the gradient of the magnetic field which is in 
the direction of B. It should be remembered that this same force is 
exerted on a particle leaving a region of more intense magnetic field, 
in which case it is then an accelerating force. 

These particles will be subjected to drifting motions due to external 
forces exerted by electric fields, the gravitational field, and inhomo- 
geneities in the geomagnetic field. We will cover these very briefly. 

If the particle moves in the presence of an electric field, E?, with 
a component perpendicular to B* then a new velocity, W*, is defined such 
that 


— > 

v 




= W + 


E X B 
B 2 ’ 


Since E and B are assumed constant. 


m 


dW 

dt 


q W X B. 


The drift velocity then results from 
the drift velocity, Wj, may be found 


the 

by 


— > 

E X B 

— term. 

B 2 


In other words. 


W d = 



In the same respect, the gravitational force on a particle having a 
component mg^ perpendicular to B yields a drift velocity 
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W d = 


The variations in the magnetic field which the particle encounters will 
also cause a third drift velocity which takes the form. 


»d - <!4 «1 + «n) 


where R is the radius of curvature of the geomagnetic line of force. 


III. MAGNETOSPHERE 

Since particles are trapped in the geomagnetic field, perhaps some 
elaboration on the earth's magnetosphere is necessary. The earth's 
magnetic field forms a definite "magnetic pocket" around us which has 
become known as the magnetosphere. The boundary of the magnetosphere 
is determined by the solar wind which "blows" radially outward from the 
sun. At a certain point above the earth, the geomagnetic field strength 
equals the kinetic energy of the solar wind. In that region, there is 
a turbulent, chaotic, intermingling of magnetic lines of force and 
energetic particles, neither of which is dominant. This region is known 
as the magnetos pheric boundary or transition zone. In the magnetosphere, 
the geomagnetic field controls, while outside the magnetosphere, the 
solar wind is the dominant energy mechanism. 

The inner boundary of the transition zone, which is called the 
magnetopause, occurs at about 10 R e on the sunward side of the earth, 
while the outer boundary of the transition zone is in the form of a 
collisionless shock wave at about 3-4 R e beyond the magnetopause. 

Figure 1 illustrates the field strength F, as well as its spatial 
coordinates, 0, and 0 (also see Reference 34). As would be expected, 
the shock wave outer boundary is most clearly seen by a smoothing out 
of F, with the direction components still variable, though less so. 

The magnetopause extends around the earth in the apparent form of 
an elongated teardrop with the tail pointing directly away from the sun. 
The shock wave, however, appears as a shroud and shows no tendency toward 
convergence. Figure 2 illustrates the recent interpretation of the IMP-I 
data. 
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-THEO. IMP-1 ORBIT 1 

NOV. 27, 1963. 
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The neutral sheet extending toward the magnetic "tail" is an area 
in which a magnetic neutral surface exists. This region appears to 
contain an excess of particles so that it was previously thought to be 
a third radiation zone. At present, however, the feeling is that these 
excess particles are necessary so that their combined magnetic moments 
may balance the geomagnetic field. It has been suggested that this 
region could be the source of the high energy auroral particles. 
Perturbations occurring in the field would produce a plasma pinch effect 
which could accelerate the low energy particles in the neutral sheet to 
very high velocities. 

The configuration of "earth, radiation belt, and magnetosphere" 
is very much like that of a typical comet with the earth representing 
the nucleus, the Van Allen belts representing the coma, and the 
magnetosphere representing the comet’s tail. 

The IMP-I data show no indication of convergence of the magnetospher ic 
boundary even though the satellite apogee extended out to 200,000 kilo- 
meters. Thus, it is presumed that the magnetospheric tail extends out 
well beyond the orbit of the moon. 


IV. THE B-L COORDINATE SYSTEM 

Perhaps the best way to explain the B-L coordinate system for trapped 
ionized particles is to take a naive approach. One normally thinks of a 
magnetic field in terms of the so-called "lines of force." Each "line 
of force" is an indication of the direction of the magnetic field along 
it, and the field strength is determined by the compactness of these 
lines of force. The following sketch (Figure 3) shows this concept for 
a simple dipole magnetic field. 



FIGURE 3 


8 



If this dipole is assumed to be enclosed in a sphere, then the 
earth's magnetic field could be roughly approximated in this way. If 
points A and C are two points along the same line of force, which may 
be designated as L, the magnetic field strength, B, is much more intense 
at C than at A. If a coordinate system is constructed based upon the 
field strength and lines of force, the resulting system could be used 
to specify equal points in any field. 

B is relatively easy to specify. It is simply the field strength 
at any point in question. But how is L defined? For simplicity, let 
us define L as the radial distance to the line of force in question, 
measured in the magnetic equatorial plane. Thus, in Figure 3 the L 
coordinate for point A is 2R, where R is the radius of the sphere 
enclosing the magnetic dipole. The L coordinate for point C is also 
2R, and here one will notice that point C does not lie at 2R radial 
distance from the center of the sphere. This discrepancy is accounted 
for by the different value of B at point C as compared to point A. 

It is interesting to realize that many satellites in near-circular 
orbits which carry experiments to determine the Van Allen radiation belt 
environment use the term L as an indication of the geomagnetic latitude. 
As can be seen from Figure 3, this is an acceptable procedure. 

Although this has been a very crude approximation, it is suitable 
in terms of the function of the B and L coordinate system. Actually, 
though B is as defined in the previous discussion, the geomagnetic field 
is not adequately represented by a dipole field. In fact, a deviation 
of 3 per cent of the calculated field at a point from the actual field 
may result in an error of more than an order of magnitude between the 
calculated flux and the measured flux. For this reason, the earth's 
magnetic field is usually determined from a spherical harmonic expansion, 
which gives an adequate representation of the geomagnetic field. 

The term L is a function of B and the integral invariant, I. If 
we denote the "mirror points" as M n and M s for the northern and southern 
hemispheres, respectively, and ds as a segment of path length along the 
line of force connecting these "mirror points," then 



Ns 
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This forms a shell along which particles are constrained to move. This 
shell contains particles with a given energy spectrum. Particles of 
higher energy levels will have their "mirror points" at lower altitudes 
within this shell, assuming equal pitch angles. The pitch angle, a, 
may be defined by 


a = arc cos — , 


where v t is the particle's total velocity and v„ is that component of 
velocity parallel to the magnetic field in the geomagnetic equatorial 
plane. 


V. SOURCES OF CHARGED PARTICLES 

One of the main problems which remain to be solved in explaining the 
Van Allen radiation belts is the question of the source of these high 
energy particles. The two main sources of particles apparently are the 
decay of neutrons caused by cosmic ray interactions with the atmosphere, 
and the "leaking in" of solar particles at the poles and diffusion 
through the magnetospheric boundary. The neutron source has been well 
established as a contributor of particles in the inner Van Allen belt. 

Cosmic rays are incident at the earth at about a rate of 2.5 particles 
(cm 2 sec)" 1 . These particles range in energies from a few Mev to about 10 9 
Bev. In the process of stopping these high energy particles, an average 
of about seven secondary neutrons are produced per cosmic ray particle. 
These neutrons then decay into protons and electrons which are in turn 
captured by the geomagnetic field. 

Solar particles are suspected of entering the magnetosphere at the 
polar regions, where the lines of force offer the least resistance. When 
these particles enter as a high energy plasma from a solar flare, the ring 
current is strongly enhanced. The precise way in which these particles 
are dispersed throughout the magnetosphere once they have entered is a 
subject of controversy. 

There is also the possibility that solar particles diffuse into the 
magnetosphere. Since there is always a solar plasma (in the form of the 
solar wind) surrounding the magnetosphere, particles could enter the 
radiation zones in this manner. The problem arises wherl one tries to 
explain how these particles with energies of a few hundred electron volts 
(in the solar wind) are accelerated into the Mev energy range. One might 
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suspect that the acceleration of particles could take place if there 
were a rapid change in B, say during magnetic storms. A large variation 
in B would produce a pinch effect which in turn could accelerate particles 
to high velocities. There is also a theory that particles may be accel- 
erated by an electromagnetic field exerting a force on the charged particles 
(Reference 39). 


VI. THE VAN ALLEN BELT MODEL 

The Van Allen belts consist mainly of protons and electrons trapped 
in the magnetic field of the earth. In the present report, we ignore the 
electron component of the Van Allen belts and deal only with the proton 
belt. This belt extends from about 1.1 R e out to about 8 Re with peak 
flux at ~ 1.5 Re measured at the geomagnetic equator. The proton component 
according to flux and energy spectra of the particles is described. 

The question of what should be used as a threshold value for the 
energy spectrum is one which is not easily answered. By moving the 
energy threshold up and down, one is able to behold a uniquely changing 
picture of the radiation zones in space. This is due to the fact that 
regions of high flux, low energy particles will be overlooked if the 
energy threshold is too high, whereas regions of high energy, low flux 
particles might be obscured by setting the energy threshold too low. 

Ibe energy spectrum is a complicated function when the energy threshold 
is low also. Thus, for manned space flight, we usually consider only 
protons with E > 40 Mev. For unmanned satellites, the problem is even 
more complicated. Since many experiments and supporting equipment may 
be exposed directly to the flux of particles incident upon them, it is 
necessary that we have a reliable radiation environment for consideration 
in designing the vehicles. Obviously, then, we need to study the entire 
energy spectrum of trapped particles. 

This problem would be challenging enough even if the trapped particles 
were temporally stable. However, this is far from true, especially in the 
case of outer zone electrons. These particles undergo spatial as well as 
temporal variations in flux by as much as an order of magnitude in a matter 
of a few hours, and for the coup de grace, even the energy spectrum at a 
particular point in space has been observed to change radically in a 
short time. 

It appears then that, if we are to construct a model of low energy, 
high flux particles in the radiation zones, it should be very simple. 

Then as our knowledge of the causative effects increases, we may update 
our model accordingly. 

It would seem extremely inadvisable at this time to attempt to construct 
a complex model based upon our present concepts of the Van Allen radiation 
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zones and their variations. Some of the problems which must be answered 
first are (1) how do trapped particles react to changes in the earth 1 s 
ring current, to changes in solar activity, and as the sun-earth-particle 
angle varies longitudinally, (2) how does the earth 1 s magnetic field react 
to these same factors, and (3) to what other, if any, perturbing forces 
are trapped particles subject? 

Let us now elaborate somewhat on possible preliminary models of the 
Van Allen radiation belt. 

One of the primary experimental problems in measuring the particles 
incident upon a counter in space is that of defining the type of particle 
which causes a particular counter to count. For instance, some frequently 
used counters on board experimental satellites are listed below (Table I) 
with the particular types of radiation to which they are sensitive. 


TABLE I 

Type of Detector Shielding Particles to which Sensitive 


Anton GM 302 

265 mg/cm 2 Mg 

Protons 

(E > 

23 Mev) 


4- 400 mg/cm 2 steel 

Electrons 

(E > 

1.6 Mev) 

Anton 213 

1.2 mg/ cm 2 mica 

Protons 

(E > 

500 Kev) 

(Directional) 


Electrons 

(E > 

40 Kev) 


48 mg/ cm 2 A1 

Protons 

(E > 

4.5 Mev) 



Electrons 

(E > 

230 Kev) 


1 . 2 mg/ cm 2 mica 

Protons 

(E > 

500 Kev) 


plus sweeping 
magnet 

Electrons 

(E > 

200 Kev) 


If a large array of detectors could be carried on board a satellite 
with a highly eccentric orbit and a reasonably long lifetime, a broad 
picture of the Van Allen belts could be obtained. 
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The Van Allen inner zone is composed of protons with peak 
intensity at about 3600 kilometers above the geomagnetic equator. The 
inner zone is fairly stable. This is to say that the flux of particles, 
as well as the energy spectrum, is constant over a long time period. 

However, one might suspect that eventually a long-term variation will 
show up in the inner zone, resulting from two effects. First, the cosmic 
ray flux incident upon the earth's upper atmosphere is about a factor of 
two higher at minimum solar activity. Therefore, more albedo neutrons 
will be available to the magnetic field for decay and subsequent capture 
of charged particles (see Section V). In addition, at least in the lower 
part of the inner radiation belt, the atmosphere may be less dense 
(Ref. 25) during minimum solar activity, thus giving the particles longer 
lifetimes. The lifetimes of trapped particles depend upon their mean 
free paths, which in turn are inversely dependent upon the average density 
of the medium through which the particles pass. 

A good method of measuring the extent of the variation in atmospheric 
density at high altitudes might be to monitor the change in the decay rate 
of inner zone particles. If the decay rate increases drastically as solar 
activity increases, the increase could probably be attributed to an increase 
of atmospheric density over the particle's length. This speculation is 
based upon the change in decay rates of artificial radiation zones r.rtsd 
by Hess (Ref. 25). 

The inner zone protons were at one time thought to be fairly localized 
below about 2 R e ; however, it now appears that low energy protons will be 
found all the way out to the magnetospheric boundary on the sunward side 
of the earth. 

Figure 4 (curve M) shows the flux versus L values of Mcllwain (Ref. 1) 
for particles of energies greater than 40 Mev. The hump at 2.2 R e is due 
to the artificial radiation particles, and the area from 3 R e out results 
from estimates of upper limits of protons. Curve R is a similar, but 
hypothetical, curve constructed on the basis of experimental data gleaned 
from several sources (Ref. 1 , 3, 6, 7 , 8, 11 , 13, 35). The second curve, 
however, assumes a low energy threshold level of 500 Kev for protons. 

Curve R does not take into account the secondary proton peak shown on 
curve M, and includes an exponential interpolation of flux With height. 

Again, it should be pointed out that in the construction of these 
graphs no attempt was made to account for temporal variations. In the 
same sense, the spurious anomalous regions will not be considered. 
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L (Earth Radii) 



FIG. 4. FLUX OF PROTONS IN THE 
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GEOMAGNETIC EQUATORIAL PLANE 


Log |Q Flux [n (> E)] 



FIG. 5. ENERGY - FLUX SPECTRUM FOR PROTONS AT VARIOUS L VALUES 
IN THE GEOMAGNETIC EQUATORIAL PLANE 
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Figure 5 illustrates the flux-energy curve we have assigned to these 
flux curves found in Figure 4. A family of curves was constructed which 
seemed to be best indicated by the various data sources mentioned in 
conjunction with the construction of Figure 4. The energy spectrum has 
been found to vary spatially as well as temporally, but no attempt is 
made here to account for this variation. 

A least squares polynomial fit to these curves yields a method of 
finding the precise value of N (> E) for any given E and L value specified. 
The data used here are again restricted to the magnetic equatorial plane. 
This method, programmed by Miss Linda Flenker, is found in Appendix A. 

We see from our previous discussion that the proton belt is subject 
to the geomagnetic field. At the nadir, this field is much changed from 
the sunward side of the earth. Thus, the proton belt will not extend 
to the magnetopause. Instead, there will be a sharp cutoff of protons 
beyond about 5 R e . This is due to the distortions of the field by the 
elongation of the magnetosphere through the action of the solar wind 
(see Section III). Because of this distortion, the particles will not 
move along lines of force defined by the integral invariant, but instead 
will be subject to perturbations caused by the stretched lines. 

In this region of the magnetosphere, the field seems to react strongly 
to such things as solar flares, enhanced ring currents, etc. The effect 
is that these perturbations serve to keep the magnetosphere beyond about 
5-8 R e relatively clear of trapped protons on the side of the earth away 
from the sun. 

Electrons in this region do appear but their occurrence is, at the 
present time, unpredictable. The electron component of the radiation 
belt will be the subject of a later report. 
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APPENDIX A 


METHOD OF CURVE FITTING BY LEAST SQUARES POLYNOMIALS 

This method gives a single equation representation of a family of 
curves (i.e., a means of obtaining a function in terms of more than one 
independent variable). 

Given: E (Energy in Mev) ; L; N (Flux) 

To Find: N = f (E, L) 

First, 


L = 1.5, 

In 

Ni 

= 

H 

O 

+ 

a n E 

+ 

a X2 E 2 

+ 

a X3 e 3 

+ 

a x4 E 4 

L = 2.0, 

In 

n 2 

= 

a 20 

+ 

a 2x E 

+ 

a ^^E 2 

+ 

a 2 3 E 3 

4 

a 24 E 4 

L = 3.0, 

In 

N-, 

J 

= 

n 

'J w 

4- 

a--E 

4- 

a-J? 2 

! 

a^^E 3 

4 

33 _ 

L = 4.0, 

In 

n 4 

= 

a 40 

+ 

a 4i E 

+ 

a 42 E 2 

+ 

a 43 E3 

+ 

a 44 E 4 

L = 5.0, 

In 

n 5 

+ 

a 50 

+ 

a 51 E 

+ 

a 52 E 2 


a 53 E3 

4- 

a 54 E 4 


A fourth degree curve was used since higher degrees increased the 
accuracy only slightly. The errors incurred in the fourth degree curve 
fit are in the hundredths decimal place. 

This curve fit gives the coefficients (a Q , a x , ... a 4 ) for which 
the equation 

In N = a 0 + a x E + a^ 2 + a 3 E 3 + a 4 E 4 

is satisfied. In this first step In N is obtained as a function of E, 
because the coefficients (a Q , a x , ... a 4 ) are a function of E. 

Second , 

a 0 's = A 0 + A X L + A^L 2 + A 3 L 3 

a x ' s = B 0 + B X L + B^L 2 + B 3 L 3 

a 2 's = C 0 + C X L + CgL 2 + C 3 L 3 

a 3 * s =Dq + D x L + DgL 2 + D 3 L 3 

a 4 's = F 0 + F X L + F^L 2 + F 3 L 3 
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A third-degree polynomial was fit in this case with errors occurring 
in the hundredths decimal place. 

These final coefficients (A, B, C, D, F) (see Table II) are a function 
of both E and L; and are obtained by fitting the coefficients 
(a 1Q> a xi , ... a 54 ) horizontally with respect to L. 

Third , 

In N = (A 0 + AjL + A^L 2 + A 3 L 3 ) 

+ (B 0 + B X L + B 2L 2 + B 3 L 3 ) 

+ (C 0 + C X L + cyL 2 + C 3 L 3 ) 

+ (D 0 + DjL + DgL 2 + D 3 L 3 ) 

+ (F 0 + F X L + F^L 2 + F 3 L 3 ) = X 


A 0 = +.249547875 + 02 


B 0 = -.295686833 - 00 


C Q = +.574196303 - 02 


D 0 = -.816146665 - 04 


F 0 = +.462310907 - 06 


TABLE II 


> 

H 

11 

-.694389747 

+ 

01 

a 3 - 

-.785188961 

- 

01 

Bi = 

-.513692717 

- 

01 

B 3 = 

-.126209083 

- 

02 

Cl = 

+.13114399 

- 

02 

C 3 = 

+.288122124 


05 

= 

-.963839773 

- 

05 

D 3 = 

+.585676357 

- 

06 

Fl = 

-.255271946 


07 

F 3 = 

-.698480307 

- 

08 


A 2 = +.125157848 + 01 
B 2 = -.224835894 - 03 
C 2 = -.224835894 - 03 


D 2 = -.144865679 - 05 


F 2 = +.392723857 - 07 
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where 


In N = f (E, L) 


to obtain 


N 

H - e £ ®> L > = e X 


Input values : 

Each input value consists of an E value and its corresponding L 
value. 


xo run ttiis program tne programmer needs the equation given in 
step 3, the coefficient values for this equation (the coefficients are 
constants for all runs), and a set of E and L values. 

Output: 

On print out sheet 


Name Value 

T Corresponds to L 

S Corresponds to E 

U Corresponds to In N 

EXPU Corresponds to N 
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